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Analysis of Superconductivity in d-p Model on Basis of Perturbation Theory 
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We investigate the mass enhancement factor and the superconducting transition temperature 
in the d-p model for the high-T c cuprates. We solve the Eliashberg equation using the third- 
order perturbation theory with respect to the on-site Coulomb repulsion U. We find that when 
the energy difference between the d-level and p-level is large, the mass enhancement factor 
becomes large and T c tends to be suppressed owing to the difference of the density of state 
for d-electron at the Fermi level. From another viewpoint, when the energy difference is large, 
the d-hole number approaches unity and the electron correlation becomes strong and enhances 
the effective mass. This behavior for the electron number is the same as that for the f-electron 
number in the heavy fermion systems. The mass enhancement factor plays an essential role in 
understanding the difference in T c between the LSCO and YBCO systems. 

KEYWORDS: d-p model, high-T c cuprate, mass enhancement factor, transition temperature, third-order 
perturbation theory 



Currently, superconductivity in high-T c cuprates is be- 
ing intensively investigated. For example, from the point 
of view of strong coupling, the t-J model has been in- 
vestigated. 1 Also, Hubbard model has been investigated 
using the variational Monte Carlo method 2 and the quan- 
tum Monte Carlo method. 3 On the other hand, from the 
point of view of weak coupling, the Hubbard model has 
been investigated using the fluctuation-exchange approx- 
imation and the perturbation theory. In the latter case, 
the nature of superconductivity in the high-T c cuprates 
has almost been clarified on the basis of the nearly 
antifcrromagnetic Fermi-liquid theory. 4 ' 5 However, we 
have not yet explained the observed differences be- 
tween the high-T c systems, YBa2Cu307_a(YBCO) and 
I^-xSr^CuO^LSCO) and so on, in particular, the prin- 
cipal reason why the transition temperature T c observed 
in LSCO is relatively low. The nuclear quadrupole res- 
onance (NQR) experiment by Zheng et al. shows that 
when the ratio of the d-hole number (rid) to the p- 
hole number (n p ), rid/2n v , is large, T c is suppressed. 
6 Actually, LSCO has a large ratio, as compared with 
YBCO. Also, in the specific heat experiment, the 7-value 
per mole in LSCO is approximately as large as that in 
YBCO. 7 - 8 Considering that YBCO has two Cu0 2 lay- 
ers and three Cu atoms in the unit, we find that the 
7-value per layer in LSCO is large. The effect of the 
strong mass enhancement can be also seen in the nu- 
clear magnetic resonance (NMR) relaxation rate. From 
the systematic study of (TiT) -1 , we can see that the ef- 
fective Fermi energy in LSCO is rather smaller than that 
in YBCO. 9 Thus, the strong electron correlation leads to 
a large mass enhancement, and then the T c in LSCO is 
reduced. This proposal is very important to clarify a re- 
lation to the other materials, such as the heavy fermion 
systems, and to promote further progress on the unified 
picture in the strongly correlated systems. In this study, 
using the third-order perturbation theory, we investigate 
the mass enhancement factor and T c in the d-p model 
for high-T c cuprates, and clarify the difference in T c be- 
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Fig. 1. Lattice structure of d-p model, i^p and t pp are the hop- 
ping integrals. 



tween LSCO and YBCO. Husimi discussed it in the same 
framework and concluded that it originates from the dif- 
ference in the effective on-site Coulomb repulsion U. 10 
We would like to stress here that the mass enhancement 
factor plays an essential role. 

Since the Cu0 2 plane is essential for superconductivity 
in high-X^ cuprates, we can consider the lattice structure 
with only the Cu0 2 network shown in Fig. 1. This struc- 
ture has the d x 2_ y 2 orbitals on Cu sites and the p x and 
p y orbitals on O sites in the primitive cell. For simplicity, 
we consider only the hopping integrals td p and t pp shown 
in Fig. 1. This is the so-called d-p model. In this study, 
we investigate in detail the nature of superconductivity 
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Fig. 2. The unbroken and broken lines arc the Fermi surfaces in 
the case of T = 0.01, i pp = 0.30 and n = 4.90 for e d - e p = 1.0 
and — £ p = 3.0, respectively. 
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in the d-p model. The model Hamiltonian is written as 



H — Hn + H\, 



(1) 



where, 



Here, e,i and e p include the chemical potential /i. The 
essential parameter is the level splitting Ed — e p . The off- 
diagonal parts are given by 
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-2t dp sin-^, (i = x,y), 



C p k = U pp sm — sm^-. 



(2) 



The second term denotes the on-site Coulomb repulsion 
between the d-electrons, 

U -J2J2 d l+ q A'- q l d k'l d kV (3) 
k q#0 



N 



We set td p = 1.0 as an energy unit, and we also fix 
t pp = 0.30 and n = 4.90 to reproduce the Fermi surfaces. 
In Fig. 2, we show the Fermi surfaces for Ed — £p — 1.0 
and Ed — £ P = 3.0 in the case of T — 0.01. We see that 
the Fermi surfaces coincide with those of the high-T" c 
cuprates, YBCO and LSCO, respectively. In this model, 
the electron number n < 5.0 represents the hole-doped 
case, and n > 5.0 represents the electron-doped case. We 
fix T = 0.01 except when we discuss T c . Figs. 3 and 4 
show the band structures and the density of states for 
d-electrons pd(s) for Ed — £ P = 1.0 and Ed — £ P = 3.0, 
respectively. We see that when Ed — e p is large, pd(s = 0) 
at the Fermi level is large. 

The bare Green's function &°\k) is written as 



G^(k) = (iu } J-H )- 1 . 
We define the matrix elements of &°\k) as 



(4) 



«L(*) 

,(0) 



GZ(k) G ( Z,(k) G^k) 
to) 



V G£' d (A) G^ Py (fc) 



(0) 



We need G^J (k) to only describe the normal self-energy 



Fig. 3. Band structure and density of states for d-electron p<j(e) 
in case of — e v = 1.0. 



2 

1 


-1 



c5 -2 
c 

CD 

-3 
-4 
-5 
-6 









I I I I I 

DOS 
























i i i i i 



X 



0.2 0.4 0.6 0.8 1 

r Pd 



Fig. 4. Band structure and density of states for d-clectron Pd(e) 
in case of — e v = 3.0. 



S(fc) and the anomalous self-energy A(fe), since the in- 
teraction exists only between the d-electrons. 

Sdd(fc) 
S(jfe) = ( 




The Green's function is given by 

G(k) = (iuj n l-H -t(k))- 1 . 



(5) 



In this study, the chemical potential /j, is determined so 
as to fix the total electron number, n. 



n = 2 §E Tr ^( fc ) = 2 §E Tr ^)- (6) 

k k 

We apply the third-order perturbation theory with re- 
spect to U. The normal self-energy is given by 

Sdd(fc) = ^[C/ 2 xo(fc-fc')G d °d ) (fc') 

k' (7) 

+ U 3 ( X 2 (k-k') + ^(k + k'))G^(k')}. 
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Fig. 5. £d — £p dependence of nd/2n p . 

Since the first-order term is constant, we can include the 
first-order term in Ed — e p . Here, 

xo(?) = -^E G d°d ) w^ d ) (?+fc), 

(8) 

^w = -§E G d d ) (fc)4°d ) (?-^)- 

k 

We also expand the effective pairing interaction up to 
the third-order terms with respect to U. 
For the spin-singlet state, the effective pairing interaction 
is given by 

V(k; k') = V RPA (k; k') + V V ertex(fc; fc')> ( 9 ) 

where 

V RPA (k; k') = U+ U 2 Xo (k - k') + 2U 3 xl(k - k'), (10) 
and 

Wertex(fc; k 1 ) = 2(T/N)Re [ ^ 4d (*i) 

hi 

x ( Xo (fc + fci)-^o(fc + A :i ))G| 1 d ) (fc + fci-fc , )t /3 • 

(11) 

Here, Vrpa(&, is called the RPA terms and 
Vvertex(& ; is called the vertex corrections. Near the 
transition point, the anomalous self-energy A(fc) satis- 
fies the linearized Eliashbcrg equation, 

A max A(fc) = E V ^ k')\G dd (k')\ 2 A(k'), (12) 
k> 

where A max is the largest positive eigenvalue. In this 
equation, the temperature with A max = 1 corresponds 
to T c . The symmetry of the superconductivity obtained 
here is the spin-singlet d x 2_ y 2 wave. 

We take 64 x 64 fc-meshes for the first Brillouin zone 
and 2048 Matsubara frequencies in the numerical calcula- 
tion. First, we show the physical quantities in the unper- 
turbed state. In Fig. 5, we show the Ed — e p dependence 
of nd/2n p . When Ed — £ p is large, nd/2n p is large. There- 
fore, in the d-p model, the difference between LSCO and 
YBCO is represented as the difference of Ed — e p . That is 
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Fig. 6. Ed — £p dependence of pd(0). 
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Fig. 7. Ed ~ E p dependence of n^. 



N 




Fig. 8. Ed — E p dependence of z^ 1 = {z^ 1 (k)). 

to say, Ed — £ P is large for LSCO and small for YBCO, in 
this model. In Figs. 6 and 7, we show the Ed — £ p depen- 
dences of pd{0) and na- When £d — £ P is large, these quan- 
tites also become large, and rid approaches unity. This 
indicates that LSCO is located in the relatively strongly 
correlated regime as compared with YBCO. 

Now, let us evaluate the mass enhancement factor us- 
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Fig. 9. £d — £p dependence of T c . 



ing the third-order perturbation theory. 

Fig. 8 shows the — e P dependence of the average 
z^ 1 = (z~ 1 (k)). Here, (• • •) represents the average over 
the momentum space. With increasing Sd — e p , the mass 
enhancement factor increases. This is consistent with the 
finding that the mass enhancement of LSCO seems to 
be larger than that of YBCO, as mentioned above. 7,8 
This is due to the difference in Pd(0). From another view- 
point, when the value of Ed — £ p is large, rid approaches 
unity, and the electron mass is strongly enhanced by the 
strong electron correlation, similarly to the Mott transi- 
tion. Such behavior has been markedly observed in the 
Ce-based heavy fermion systems. 

Next, we discuss the superconducting transition tem- 
perature T c . From the experimental results, for large 
£d — £p, T c is relatively low. In Fig. 9, we show the Sd — £ P 
dependence of T c . T c for U = 4.0 increases as a function 
of Ed — e p , and for U = 6.0 is almost unchanged. Al- 
though the normal self-energy correction effect markedly 
suppresses T c , we cannot derive the results observed in 
the experiments within the present calculation. On the 
basis of the following reason, however, we can expect 
that by calculating the higher-order terms, we will ob- 
tain a better dependence of T c on Ed — e p . In Fig. 9, the 
filled (open) circles and squares denote T c for U = 4.0 
(U = 6.0) with and without the normal self-energy cor- 
rection, respectively. The difference is mainly due to the 
effect of the mass renormalization. The larger Ed — £ P 
is, the more this effect increases. This tendency is con- 
sidered to be more marked in the higher-order terms of 
the normal self-energy. On the other hand, higher-order 
terms of the pairing interaction were evaluated by No- 
mura et al. in the Hubbard model. The results show 
that the convergency for the pairing interaction is very 
good for the spin-singlet pairing near the half-filling. 11 
Namely, the contribution of the fourth-order terms to 



the pairing interaction is small compared with those of 
second-order and third-order terms. Here, we can expect 
the same trend. Thus, with the inclusion of the fourth- 
order terms in the Fliashberg equation, T c will become 
relatively low for large Ed — e p . This prediction is an im- 
portant problem to be confirmed in the future. 

In conclusion, we have investigated the mass enhance- 
ment factor and the superconducting transition temper- 
ature in the d-p model for the high-T c cuprates. We have 
solved the Fliashberg equation using the third-order per- 
turbation theory with respect to the on-site repulsion U. 
We find that when £d~£ P is large, the mass enhancement 
factor becomes large and T c tends to be suppressed ow- 
ing to the difference in pd(0). From another viewpoint, 
when the d-hole number approaches unity, the electron 
correlation between d-holes (electrons) becomes strong 
and the effective mass increases. In fact, LSCO with d- 
hole number near unity shows strong mass enhancement . 
Here, we consider that T c is given by the renormalization 
factor z and T* as 

T c * zT c *, (14) 

where T* is the critical temperature determined by 
the calculation without any renormalization due to the 
normal self-energy correction. LSCO possesses small z 
and exhibits low T c . Also, in the heavy fermion systems, 
when the number of f-electrons approaches unity, the 
effective mass is large and T c is suppressed. This is 
the important unified theory which holds for all the 
strongly correlated electron systems from cuprates to 
heavy fermions. Thus, our calculation shows that, in 
order to systematically discuss the physical quantities 
such as T c and the mass enhancement factor, we need to 
use the d-p model. The Hubbard Hamiltonian is insuffi- 
cient to represent the difference among material systems. 
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